Exact renormalization group approach to a nonlinear diffusion equation 
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The exact renormalization group is applied to a nonlinear diffusion equation with a discontinu- 
ous diffusion coefficient. The generating functional of the solution for the initial- value problem of 
nonlinear diffusion equations is first introduced, and next a new regularization scheme is presented. 
It is shown that the renormalization of an action functional in the generating functional leads to an 
anomalous diffusion exponent in full order of the perturbation series with respect to a nonlinearity. 
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The renormalization group (RG) is a powerful tool to 
reveal universal behavior of various systems, including 
quantum field theories and statistical mechanics Its 
basic idea lies in the coarse-graining of short-distance de- 
grees of freedom, which causes redefinition of parame- 
ters governing the long-distance physics of the systems 
under investigation. In spite of its conceptual simplic- 
ity, there exist in the RG methods many techniques to 
solve with, and numerous applications have been made to 
equilibrium or near-equilibrium systems 0, Q . In partic- 
ular, exact RG (ERG) (also called nonperturbative RG 
or functional RG) techniques and some approximations 
based on them P, 0, S S Q have been attracting much 
renewed interest to reveal nonperturbative phenomena in 
field theories la ^statistical mechanics 

9]. 

and condensed 

matter physics 0, 0, . 

On the other hand, Goldenfeld et al. [l3L fl4| have ex- 
tended the RG methods to systems far from equilibrium. 
They have demonstrated that there exists a deep rela- 
tionship between the RG and the intermediate asymp- 
totics method [l5| in the study of the nonlinear partial- 
differential equations for nonequilibrium systems. Their 
idea has attracted much interest, and the RG approach to 
nonlinear differential equations has been developed [lfij . 

In this paper, we apply the ERG method to a nonlinear 
diffusion equation called Barenblatt equation [l^. This 
equation has a discontinuous diffusion coefficient; this 
discontinuity makes perturbative expansion more com- 
plicated if one proceeds to higher order computations. 
We show in this paper that Polchinski equation, a ver- 
sion of the ERG equation, is very efficient even for such 
a nonlinear diffusion equation. It turns out that we can 
indeed solve the equation for all order in the perturba- 
tion series. The solution leads us to the full anomalous 
diffusion exponent. 

Let us start with the following nonlinear diffusion equa- 
tion called Barenblatt equation; 



d t u{x, t) - D{u)d 2 x u(x, t) = 



(1) 



with an initial condition u(x,t = 0) = q5(x), where 
D(u) = k[1 + g6(— adtu)] denotes a nonlinear diffusion 
coefficient with a being a positive constant which makes 
adtu dimensionless (below, we set a = 1, for simplicity). 



Here, 6(x) = (1) for x < (x > 0) stands for the 
step function. The dimensionless constant g controls the 
nonlinearity of the diffusion coefficient. The Barenblatt 
equation describes the filtration of a compressible fluid 
through a compressible porous medium which can be ir- 
reversibly deformed. Goldenfeld et al. 0] have obtained 
asymptotic behavior of the solution by solving this equa- 
tion via an iteration scheme corresponding to the pertur- 
bative expansion with respect to g. This perturbation 
gives rise to divergences: Their basic idea is introduc- 
ing a renormalization scheme which render the solution 
finite and deriving an anomalous diffusion exponent as 
an anomalous dimension in the RG language. Though 
they have successfully obtained the leading correction of 
the diffusion exponent, their method seems difficult to 
extend to higher order due to the discontinuous step- 
function nonlinearity. We will present a new renormal- 
ization scheme for the initial-value problems of nonlinear 
diffusion equations. 

To be specific, we introduce a generating functional of 
the solution for the Barenblatt equation and a modified 
propagator with a short-time cutoff to render the solu- 
tion finite. To this end, notice 0,0] that the solution of 
Eq. can be written as u(x,t) — J T>(f><fi(x,t)5(dt(j> — 
D((j))d 2 (j))S((j)(x,Q) — u(x, G)). This expression can be 
rewritten by a functional integral if the derivative with 
respect to t is interpreted as a forward difference op- 
erator. Namely, using the Fourier transformation for 
the delta function, we reach u(x,t) = (<f>(x,t)) = 
i J Vcj/Dcf) (/>(£, i)e lS with S being an action functional 



s 



dt 



dx 



<pA~ V - g(/>6(-dt<t>)Kd%.<t> - <pJ ,(2) 



where A(x, t) = e~ x /( 4Kt ) jyfAmrti denotes the diffusion 
propagator and the generating functional Z is defined 
by Z = fT><f)D(f> e lS , as usual. The field J in the last 
term is defined by J(x,t) = u(x,0)8(t) which controls 
the initial value of u. In what follows, we examine the 
case u(x,0) = qS(x), as mentioned below Eq. QJ, but it 
should be stressed that the generic initial-value problem 
can be treated similarly. 

As discussed by Goldenfeld et al. perturbative 
calculation diverges with the initial condition specified 



2 



above. To regularize the solution, they have introduced 
an initial distribution with a finite width such as u(x, t = 
0) = er x ^ 2l ' /V2irl 2 . We instead introduce the short- 
time cutoff e for the propagator to formulate the ERG for 
the present nonlinear diffusion equation. To be specific, 
we define a modified propagator as 



where S £ is defined by 



-iS e [J-A > +<t> < ,A>-J+<f><] 
= p-iSe [i&/SJ+J-A > +4> < ,i6/5J+A > -J+<p < ] 



(8) 



A e (x,t) = 6(t - s)A(x,t). 



(3) 



One can easily check that this propagator indeed gives 
a finite solution in the perturbation theory, applying it 
to the calculations by Goldenfeld et al. |lj|. This regu- 
larization scheme can be used not only in the Barenblatt 
equation but also in generic diffusion problems. 

Having defined the generating functional and the mod- 
ified propagator, we next derive an ERG equation for the 
action functional Eq. (0 ■ Using the propagator @ with 
a cutoff £o and introducing a source term, we start with 
the generating functional, 



This equation implies that if we expand the exponential 
in r.h.s and make all the derivatives 5/5 J and 5/5 J in 
S £o act on J and J in the right S £a , we reach some S e as 
a functional of $ = J • A> + c/) < and $ = A> • J + <j) < . 
Instead of carrying out such calculations, however, we 
can alternatively determine the functional S e by noting 
that Z £ obeys 



dZ £ 
de 



5J 



de 



5J 



(9) 



which follows from Eq. 10. Substituting Eq. Q into 
Eq. ©, we obtain the following Polchinski RG equation, 



Z[J J] = 2?02?^ e i ^' A £o 1 '^ _iSe o^'^~ iJ '^ ,_i *' J (4) _ 5S £ ^ dA > ^ 5S £ 

de 



5S e dA> 
c><f> de 



5S e 
5$ 



itr- 



dA- 



5 2 S £ 



de <5<M$ 



(10) 



where the bare action of the nonlinear term is S E0 [<fi, <fi] = 
g<f> ■ 9(—dt(f)Kd 2 (t). At the end of the calculations, we 
must set J(x,t) — q5(x)5(t) to obtain the solution for 
the initial-value problem of the present equation. Here, 
the symbol a ■ b implies a-b = J Q dt J_ dxa(x, t)b(x, t). 
Next, we introduce a new cutoff e (> eo) and divide the 
propagator into two parts A eo = A> + A< , where 



Next task is to determine the functional S s by solving 
this equation. To this end, let us first consider the $ 
dependence of the functional S £ . The bare S £o contains 
only the first order term in but in the process of the 
renormalization, Eq. (|l()(l yields the zeroth order term in 
S e - To be concrete, let us denote 5 e [$, $] = $ • H e [$] + 
F[$\. Substituting this into Eq. l(TU|) . we find 



A> = [0(t - so) - 6{t - a)] A, 
A< = 8{t-e)A. 



(5) 



Here, > and < imply the short-time and long-time 
modes, respectively. Separating^also 6 and <f> into two 
fields 4> = 4*> + < \>< an d 4> = 4>> + 4>< enables us to 
rewrite the generating functional as 



de 
dF £ 

de 



5H £ dA~ 



5$ de 
5F £ dA> 

(5$ de 



itr- 



dA > 
de 



(11) 

SH £ , s 

ii (12) 



Z[J,J] cx V^V^e 



i^K-A- 1 - 



'<Z S [J,J; 0<,0<], 



Z £ [J, J; = / VfeVfee^ 



■id> > ■ A^ 



x e 



-iSe [<t>>+4><,(t>>+<t><\-iJ-(4>>+4><)-i(4>>+4><)-J 



(0) 



up to a proportionality constant. The field </>>, (j) > and 
(j) < , <j) < can be identified as fields describing short-time 
and long-time modes, respectively. Integrating out the 
short-time fields, we will next derive an effective action 
describing long-time modes. 

^ Changing^ the integration variables (j) > and 0> into 
4> > = 4> ~jt>< an d 4>> = 4> — 4><i an d integrating over 
the fields <j> and cj> in Eq. © yields 

2 £ [J, J; <P <,(/)<} 

q}^q}^ e Ki-i<)-A> 1 -{4>-4><)-iS S0 [i,cp}-iJ-ct>-ii-J 

-iJ ■ A> • J —i<j> < ■ J — ij ■ 0< — iS c [J- A> + 0< ,A> ■ J+<ji < ] 



with the bare functions -f^J^] = g0(—d t ^)ndl^ and 
F £o [<I>] = 0. The RG equation for the -He-term, which 
determines the solution of the Barenblatt equation, is 
closed. Furthermore, it has no loop corrections. Nev- 
ertheless, the initial-value problems are still nontrivial, 
since after obtaining H £ , we must set J = q5(x)5(t) and 
determine the e-dependence. 

To solve the functional equation (|TT|l . we assume the 
form of H £ as 

/>oo />oo 

H e [$](x,t) = ds dyV £ [dMx-y,t- S )} K d 2 ^(y,s) 

JO J -oo 

= V e [$\'K$"(x,t) (13) 

with a certain unknown function V £ , where we have de- 
noted <9 t $ = $ and d 2 <& — $" for simplicity. Substituting 
this into Eq. we have 



d £ V £ [^} = V £ [$] ■ d £ KAl ■ V e [$\. 



(14) 



The bare function is given by Ve [$](a; — y,t — s) = 
—g9(—<&(x,t))5(x — y)5(t — s). This equation can be 
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solved if V e is expanded in power series of g such that with 5(f)(\ / 2Ktuj, t) = (f> < (\/2Ktuj, t) — 4> < (0, 0) and f{ui) 



V, = El 

obeys 



1 g n V e (n) . Actually, we find that each term 



= E 



^ (ni) [$]-d e KA'^-I/(™ 2 )[$]. (15) 



(i) 



and 



From this equation, it follows that d e V e 
hence, it turns out that Ve is not renormalized; that 
is, 14 (1) [<j>] = V[$], where 

V[$\(x - y,t - s) = 6(-$(x,t))5(x - y)8(t - s). (16) 

This enables us to calculate the higher order solutions 
by substituting Eq. 1)16(1 into Eq. (|15|l and by solving it 
successively order by order: 



V e (n) = V (-kA" • Vf 



(17) 



Thus, we have determined the functional ff e [ < i ) ] as infi- 
nite power series with respect to g. In passing, we briefly 
discuss the solution for F e . Expanding similarly the RG 
equation 1121) in power series of g, we find that the solu- 
tion of Eq. 1)12)1 is F e — because of the fact that the 
bare F ea — as well as that the second term in the r.h.s 
of Eq. 1)12)1 is zero due to the trace with respect to the 
time variable. Therefore, we end up with the renormal- 
ized action functional, 



S F = $ 



oo 
n=l 



n V[$] f -kA" • V[<$>}\ ■ k$" 



9' 



(18) 



with V defined by Eq. 1)16)1 . Thus, we have determined 
the renormalized action in full order in g. This shows the 
efficiency of the present approach. 

To obtain the asymptotic solution of the Barenblatt 
equation, we must set J(x,t) = q5(x)5(t) to specify the 
initial condition and calculate dominant parts with re- 
spect to e/e in the action obtained so far. The ac- 
tion decomposed into some sectors by the substitution 
of $ = J • A> + 4> < and $ = A> • J + <f> < . Then, we 
notice that among them dominant contributions are from 
the $(x,t) = <j) < (x,t) and $(ar, t) = qA(x,t) sectors in 
the asymptotic region e<Ct. 

Let us start with the first order which Goldcnfcld et al. 
have calculated within the perturbation in the iteration 
scheme. From Eq. 1)18)1 . one can obtain 




<fe0< (x, t)9{-qA > (x, i))«A" (x, t) 
0<(O,O) + gS reg (19) 



where the regular part of is defined as 



c(i) = 



dt 

T 



dus 5<p < (\/2KttJ, t) f{uj) 



(20) 



-1 e~ 



/2 



— n -j— -. Here, the regular term in Eq. 1)19)1 is not in- 
volved with the renormalization of the action since we can 
safely set £o — * 0, while the first term is relevant to the 
renormalization of q, the height of the initial distribution. 
Thus it turns out that at this order q is indeed renormal- 
ized, whereas others, especially g, is not renormalized. 
This feature actually holds even in the next order, as will 
be checked below. This implies that the present system 
is always at a fixed point because g is not renormalized, 
and hence, the anomalous dimension which is in general 
a scheme-dependent quantity is a physical observable in 
the present case. Considering these, we introduce the 
renormalization only to q and define the renormalized q 
as qn = qZ . Expanding the renormalization constant Z 

as Z = 1 + E n =i9 nz{n) with zin) = -7 (n) hi(e/£o), it 
turns out that the first order of Z reads 



7 



(i) _ 



1 



(21) 



from Eq. ()19)l . This indeed reproduces the result of 
Goldenfeld et al. 0. 

The renormalization of q introduced above is indeed 
enough to render the solution of the Barenblatt equation 
finite also in higher order. To verify this, let us next 
calculate the second order renormalized action. Due to 
similar arguments to the first order, the action 1)18(1 yields 



5i 2 



I 2k l 



t—e 



q dt dx dti / dxi 

J2sq J — J so J — \/2Kt\ 

x<j> < (x,t)KA"(x - xi,t-ti)K,A"(xi,ti) 
+reg. , (22) 

where reg. stands for regular parts of the renormalized 
action. Similar but a bit lengthy calculation leads to 



s< 2 > = 



Z^ 



2 

+reg. 

where Z^ = —^f( 2 > ln(e/e ) with 

7 (2) = -t—t ^i/M 



#0,0) + qZ^S { r 1 c l + qZ (2 U(0,0) 



(23) 



Jo T i J-i 
1 1 - vn^i e -(i-V^i)7[2(i-n)] 



l-n 



(24) 



This result indicates that the second order action cor- 
rectly includes the contributions from the first order 
renormalized action. Namely, the renormalized action 
with the source term q<f>(0, 0) satisfies gS^ + g 2 sf > + 
q<p(0, 0) = e z qS lcg + e z q<j){0, 0) up to the second order of 
g. Therefore, we expect in general that qn = q{e^/e) 1 
where 



7 



E^7 ( " 



(25) 
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The constant 7 thus obtained indeed gives the anoma- 
lous dimension of the solution for the present diffusion 
equation. 

To see this, notice that since J is not renormalizcd, 
we have u(x, t; q, £0) = u(x, t; qr, e), which tells that the 
solution is independent of e. Hence, the renormalized 
solution should satisfy the following RG equation; 



d_ 



d 

iqR ~dq~R ^ U<yX ' t; qR ' ^ = °' ^ 



where 7 is defined by Eq. Q25[l. or alternatively by 
7 = — dlnqji/dlne. On the other hand, the dimensional 
analysis requires u(x, t; q, Eq) = q/yKt$ (x/y~Kb, £oA) • 
Therefore, combining these observations, we can as- 
sume u(x,t;qn,e) = qn/y/Kt^ (x/^J~Kt, e/i) . Substitut- 
ing this into Eq. I|26() . it turns out that the relation 
$ (x/Vfii, e/t) = (e/t)' 1 1> (x/^/kF) holds. Hence, the 
asymptotic behavior of the solution for the Barenblatt 
equation is indeed given by u ~ l/t 1 / 2+7 . 

So far we have derived the anomalous diffusion expo- 
nent up to second order with respect to g. The exponent 
of first order is the same as that obtained by Golden- 
feld et at, whereas the exponent of second order l|24|) is 
estimated as 7^) ~ —0.063546 by numerical integration. 
Though this value is different from the result in Ref . , 
it coincides with that obtained later by Cole and Wagner 
given by a different integral formula derived via differ- 
ent method [hsj . Since the anomalous diffusion exponent 
should be scheme-independent in the present case, we 
believe that our result as well as Cole and Wagner's is 
correct. 

The efficiency of our method lies in the fact that one 
can compute the higher order exponent in a similar way 
above. To present the exponent of nth order, we define 
a function, 



g(ui,u) 2 ,T 2 ) 



-(wi-Vfiwa) /[2(1-t 2 )] 

V27r(l-r 2 ) : 



/2vr 



(27) 



Then, Eq. iJTSJ yields 



1 71 — 1 




I n—1 

n — / n ^'/k-o 

1 = 1 Tj J - 1 3 = l 

1 - y/nui e-d-^J'/PP-n)] 



1-Tl 



1 d 2 g(u) j ,u) j+ i,T j+1 ) 



(28) 



We expect that this anomalous exponent is exact, since 
it has been derived from one formula (|18|l . among which 
the first and second order coincides with those calculated 
by different methods. For references, we numerically esti- 
mate the exponent of third order, 7^' = —0.00314. The 
exponent obtained so far seems a good convergent series. 



In summary, we have applied the ERG techniques to 
the initial- value problem of the Barenblatt equation, one 
of typical nonlinear diffusion equations. We have de- 
rived the anomalous diffusion exponent in full order with 
respect to the parameter controlling the nonlinearity. 
This implies that the ERG approach is efficient to sys- 
tems far from equilibrium described by nonlinear partial- 
differential equations as well as to field theories and sta- 
tistical mechanics. 

The present approach would be useful to other types 
of initial- value problem for nonlinear diffusion equations. 
In particular, application to the critical dynamic of non- 
linear traveling wave is of great interest. One of well 
known examples is the KPP equation |l9l| which shows 
an interesting universal behavior in the selection of the 
front velocity. This problem has been addressed by Pa- 
quette et al. [2(|, but more detailed analysis is needed 
if one wants to understand, e.g., the universal logarith- 
mic corrections to the velocity in the pulled- front, from 
the RG point of view. Application of the ERG to such 
problems would be quite interesting. 
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